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An integral-differential equation is derived for the self-consistent (effective) 
field in the medium consisting of many small bodies randomly distributed in some 
region. Acoustic and electromagnetic fields are considered in such a medium. Each 
body has a characteristic dimension a <C A, where A is the wavelength in the free 
space. The minimal distance d between any of the two bodies satisfies the condition 
d ^ a, but it may also satisfy the condition d <C A in acoustic scattering. In 
electromagnetic scattering our assumptions are a <C A and A <C d. Using Ramm's 
theory of wave scattering by small bodies of arbitrary shapes, the author derives 
an integral-differential equation for the self-consistent acoustic or electromagnetic 
fields in the above medium. 

1 Introduction 

A general method to derive equations for the self-consistent (effective) field in a medium 
consisting of many small particles is proposed. The method is illustrated by the derivation 
of such equations for acoustic and electromagnetic waves. The novel points in this paper 
include: 

1) the method for the derivation of the equations for the self-consistent field in a 
medium consisting of many small bodies of arbitrary shapes; 

2) equations (2.29), (3.2)-(3.5), and the formulas for calculation of the polarizability 
tensors for bodies of arbitrary shapes with the desired accuracy. 

Equation (2.16) is of the type, obtained earlier in by a different argument, and in 
It is simpler than equation (2.29). This can be easily explained physically: scattering 
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Abstract 



of an acoustic wave by small, in comparison with the wavelength in the free space, 
acoustically soft body is isotropic, and the scattered field in the far-field zone is described 
by one scalar, the electrical capacitance of the perfect conductor of the same shape as 
the small body, and the scattered field is of order 0(a), where a is the characteristic size 
of the small body. If the small body is acoustically hard, that is, condition (2.20) holds 
on its boundary, then the scattering is anisotropic, the scattered field in the far-field zone 
is described by a tensor, and this field is of order k'^a?, that is, much smaller (by a factor 
of order 0((A;a)^)) than for the acoustically soft body of the same size a and of the same 
shape. Here k is the wavenumber in the free space. 

In jTj the Lorentz-Lorenz formula is derived. This formula relates the polarizability of 
the uniform dielectric to the density of the distribution of molecules and the polarizability 
of these molecules. In this theory one assumes that the polarizabilty in the matter is 
a constant vector, the molecules are modeled as identical spheres uniformly distributed 
in the space. In this case the polarizability tensor is proportional to the unit matrix, 
and the coefficient of proportionality is the cube of the radius of the small sphere times 
some constant. This, together with additional assumptions, yields a relation between the 
dielectric constant of the medium, the polarizability of the molecule, and the number 
of the molecules per unit volume. The derivation of this formula in is based on the 
equation of electrostatics. 

Our basic physical assumptions, (1.1), allow for rarefied medium, when d » A, but 
also, in acoustic wave scattering theory, for medium which is dense, when a « d « \. 
Equations (2.29) and (3.2)-(3.5), that we have derived, have an unusual feature: the 
integrand depends on the direction from y to x. This happens because of the anisotropy 
of the scattered field in the case of non-spherical homogeneous small bodies. 

A possible application of equation (2.16) is a method for finding the density of the 
distribution of small bodies from the scattering data. Namely, the function C{y) in (2.16)- 
(2.17) determines this distribution. On the other hand, this function can be determined 
from the field scattered by the region TZ. The uniqueness results and computational 
methods for solving this inverse scattering problem are developed in 0,101, and |7j (see 
also the book jHj pubhshed in 2005). 

Below we study the dynamic fields, so that the wavenumber is positive. 

Consider a random medium consisting of many small bodies Dj, 1 < j < J, J 10^^, 
located in a region TZ. Let aj be the radius of the body Dj, defined as a,- := ^ max \x — y\, 

x,y£Dj 

and a := max a,,-. We assume a <^ A, where A is the wavelength of the field in the 

i<j<J ^ 

free space (or in a homogeneous space in which the small bodies are embedded). Let 
d := min \x — y\. Assume 

x&Dj,y&Dj,l<i,j<J,i^j 

a < A, a < d. (1.1) 

We do not assume that d^ X 'm. the acoustic wave scattering, but assume this in electro- 
magnetic wave scattering. The difference in the physical assumptions between acoustic 
and electromagnetic theory is caused by the necessity to apply twice the operation Vx 
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to the expression of the type in the electromagnetic theory, where p is a vector 
independent of r and r is the distance from a small body to the observation point. 

We consider acoustic field in the above medium, and derive in Section |21 an integral- 
differential equation for the self-consistent field in this medium. The notion of the self- 
consistent field is defined in Section |21 Roughly speaking, it is the field, acting on one of 
the small bodies from all other bodies, plus the incident field. 

In Section IHl we derive a similar equation for the self-consistent electromagnetic field 
in the medium. 

Each small body may have an arbitrary shape. The key results from |4j, that we 
use, are the formulas for the 5'-matrix for acoustic and electromagnetic wave scattering 
by a single homogeneous small body of an arbitrary shape. These formulas are given in 
Section 121 and Sectional Wave propagation in random media is studied in [2]. 

2 Acoustic fields in random medium. 

Assume first that the small bodies are acoustically soft, that is m |5.= 0, where Sj is the 
boundary of Dj. The governing equation for the acoustic pressure u is 

{V^ + k'^)u = Om D' ■=R^\D, D := i Dj, (2.1) 

(2.2) 

u G ^^ (2.3) 



J 

u = on S := U Sj, 



dv \ 

ikvj^O, r := \x\, (2.4) 

where the direction u & oi the incident field is given, k = ^ = const > is the wave 
number, is the unit sphere in R'^. 
Let us look for v of the form: 



^ r ^ik\x—y\ 

vi^) = y2 9{x,s)(7j{s)ds, g{x,y):=— -. (2.5) 



^ik\x—y\ 

g{x,s)(Tj{s)ds, g{x,y):=- 

j=i 

If fjl.l|) holds, and min |a; — Sj| ^ a, Sj G Sj, then ()2.5|) can be written as 



l<i<J 

v{x) = y2g{x,Sj)Qj[l + 0{ka)], Qj '■= / o-jds, (2.6) 
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Define the self- consistent field at the point Sm G 5"^ by the formula 

J 

Ue{Sm) = Uo{Sm) + ^ giSm,Sj)Qj, 



and at any point x, such that 
by the formula: 



infj\x — Sj\ ^ a, 



Ueix) = Uo{x) + ^g{x, Sj)Qj. 

i=i 



(2.7) 

(2.8) 
(2.9) 



If (ESI) holds, then 



^g{x,Sj)Qj - g{x,Sj)Qj 

j=l j=ljjtm 



< \Ue{x)\, 



that is, removal of one small body does not change the self-consistent field in the region 
which contains no immediate neighborhood of this body. On the surface S of the body 
the total field u = 0, so the self-consistent field on 5*^ differs from u, while at a point x 
such that (ITHl) holds, u{x) = Ue{x) + o{l), if dHH) holds. 
Let us derive a formula for Qj. By (j2.2|) one gets: 



= UeiSm) + 



(Tm{s)ds 

47r|sm, - s| 



+ £r 



where 



so that 



[g{sni, s) - go{sm, s)]amds, go :-- 



4:7c\Sm - S\ 



maxl^ml = o(l) as ka — > 0. 

m 

Thus, one may neglect Em in ()2.1()|1 and consider the resulting equation 

JSm 47r|sm - s\ 



(2.10) 



(2.11) 



f2.12) 



(2.13) 



as an equation for the charge distribution on the surface Sm of a perfect conductor 
charged to the potential —Ue{sm)- Then, by ()2.6|) . 



(2.14) 
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where Cm is the electrostatic capacitance of the conductor Dm- From (|2.9|) and (j2.14j) 
one gets 

J 

Ue{x) = Uo{x) -^g{x,Sj)CjUj{sj) + o(l), (2.15) 
i=i 

as ka —>■ 0. 

Let us emphasize the physical assumptions we have used in the derivation of (2.15). 
First, the assumption ka << 1 allows one to claim that, uniformly with respect to all 
small bodies, the term o(l) in (2.15) tends to zero as ka —>■ 0. Secondly, the assumption 
d » a allows one to claim that the m— th small body is in the far zone with respect 
to the j— th body for any j ^ m. The expression under the sum in (2.15) is the field, 
scattered by j— th body and calculated at the point x, such that \x — Sj\ » a, that is, 
in the far zone from the j— th body. So, physically, the equatipons for the self-consistent 
field in the medium, derived in this paper, are valid not only for the rarified medium 
(that is, when a << X and d >> X, but also for not too dense medium, that is, when 
d » a and a « A, but, possibly, d « X. 

The limiting equation for Ue{x) is: 



Ue{x) = Uo{x) - I g{x,y)C{y)ue{y)dy, (2.16) 

where 



C{y)dy= J2 Q> (2-17) 

Sj Gdy 

and the summation is taken over all small bodies located in the volume dy around point 
y. If one assumes that the capacitances Cj are the same for all these bodies around 
point y, and are equal to c{y), then C{y) = c{y)N{y), where N{y) is the number of small 
bodies in the volume dy. 

Equation (j2.16|) is the integral equation for the self-consistent field in the medium in 
the region TZ. This field satisfies the Schrodinger equation: 

[V^ + k'^ -q{x)]ue = 0, q{x):=C{x). (2.18) 

Since Cj ~ a in (|2.17|) . and the number of the terms in the sum (j2.17p is = 0{j^), 
provided that dy is a unit cube, one concludes that ^ = 0(1), so 

d=o(a^). (2.19) 



If one had Na^ = 0(1), i.e., small bodies have nonzero limit of volume density, then the 
assumption d^ a would be violated. 

Let us now assume that the small bodies are acoustically hard, i.e., the Neumann 
boundary condition 

Oil 

— = OonS, (2.20) 
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replaces (|2.2j) . M is the outer normal to S. In this case the derivation of the equation for 
Ue{x) is more complicated, because the formula for Qj is less simple. It is proved in jl] 
that for the boundary condition ()2.20|) one gets 

Q, = ikV,P^;^^np-^ + V,Au,. (2.21) 



Here and below, one sums up over the repeated indices, Vj is the volume of Dj, Au^ is the 
Laplacean, Up := the small body Dj is located around point y, the scattered field 

is calculated at point x, /3pg is the magnetic polarizability tensor of Dj which is defined 
by the formula j3] 

Ppg:=apg{-1), (2.22) 
and a!pg(7), 7 := is the electric polarizability tensor, defined by the formula: 

Pp = apg{-f)VgeoEg. (2.23) 

Here P is the dipole moment induced on the dielectric body Dj, with the dielectric 
constant e, placed in the electrostatic field E in the homogeneous medium with the 
dielectric constant Eq. 

Analytical formulas for calculation of dpqi'y) with an arbitrary accuracy, in terms of 
the geometry of Sj, are derived in j^. These formulas are: 

m=0 ^ ' 

where 

«':=V^„.6»:=/ /^!^(£)»^,r„:=|.-*U«=^ = (2.26) 

Js, Js, \0, if p^q, 



•= / / dsdtAfp{s)Xgit) I ... I —Hti,t) . . . ^{tm-l,trn-2)dti . . . dtm-l. 




(2.26) 

(1) 



PI 

and 

PQ 



PS = (2.27) 



«P.(7) - «i?(7)l = 0(g"), 0<q<l, (2.28) 
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that is, (|2.24|) gives a convergent approximation of the tensor apg{j). 
From (Q and ^TT^i one gets 



Xp - ypduejy) 
x-y\ dyg 



Ue{x) = uo{x) + / g{x,y) 
Jn 

where V{y) and (3pq{y) are defined by the formulas 



V{y)dy, (2.29) 



V{y)dy = J2 ^,{y)V{y)dy = J] (3p,,V,, (2.30) 

where (ipqj is the magnetic polarizabihty tensor of the j — th small body, and the summa- 
tion is over all small bodies in the volume dy around point y, so that V{y) is the density 
of the distribution of the volumes of small bodies at a point y. 

Equation (2.29) is approximate, in contrast to (2.16). Indeed, if a << d then 
lima^o Ss^gdy ~ because Vj = O(a^) and the number of small bodies in the unit 
volume is O(^), so lim«^o edy ~ ^iw^y) ~ ^- However, if the Dirichlet condition 
(2.2) holds, rather than the Neumann one (2.20), then ^ = 0(1) and equation (2.16) is 
exact in the limit f — ^ with C{y) defined in (2.17). 

The novel feature of equation ()2.29|) is the dependence of the integrand in ()2.29p 
on the direction x — y. This one can understand, if one knows that the acoustic wave 
scattering by a small soft body Dj is isotropic and depends on one scalar, the electrostatic 
capacitance of the conductor Dj, while acoustic wave scattering by a small hard obstacle 
is anisotropic and depends on the tensor Ppg. 

Finally, if the third boundary condition holds: 

^ + hu = on S (2.31) 
then (see 0), if /i is not too small, one has: 

^ ^j^ij /2 32) 

where \Sj\ is the area of Sj, so that 

Ue{x) = Uo{x) - [ g{x,y)h{y)ue{y)dy, (2.33) 

where 



n 



Sj ady 



8 



3 Electromagnetic waves in random medium 



In this section our basic assumptions are: 

a « X d» X. (3.1) 

The reason for the change in the assumption compared with (1.1), where d is not neces- 
sarily greater than A, is the following: formula (3.4) below is valid if a << A and d » X, 
where d is the distance from a small body to the point of observation. This comes from 
applying twice the operation of Vx to the vector potential (see, e.g., |j4j, p. 98-99). In 
acoustic scattering the formula for the scattered field (see e.g. (2.7)) is valid if d >> a, 
while formula (3.4) (see below) is valid for d >> A and a « X. 

Let U := (^). Denote by S the 6x6 matrix which sends U into gUsc, where gUsc is 
the scattered field, g := — , r is the distance from a small body located at a point y to 
the observation point x, r = \x — y\. 

If S is known, then the equation for the self-consistent field Ue in the random medium 
situated in a region TZ, and consisting of many small bodies, satisfying conditions (3.1), 
is 

Ue=Uo+ [ g{x, y)S (y, U,{y)dy, (3.2) 



Jn 

where 



(3-3) 

as follows from the argument given for the derivation of ()2.29p . 

Let us give a formula for S^^\ assuming without loss of generality, that the origin 
is situated inside a single body Dj, which has parameters ej, fxj, and aj, (dielectric 
permittivity, magnetic permeability, and conductivity, respectively) , and drop index j in 
the formula for S^^^ . 

In j3] one can find the formulas 

EsM) = ^ (-W, [P, u']] + ./^[M, ly']) , (3.4) 
47r \eo V £o / 

Hscii^') = Escii^')] = ^ (^^W, P] + [M, u']]] , (3.5) 

where is the unit vector in the direction of the scattered wave. 
Formulas ()3.4|) . ()3.5|) can be written as 
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where the matrix S is expressed in terms of the tensors a := apq and P := f3pq, because 
P is calculated in (|2.23p . 

Mi = aij{j)V fioHj + PijV fioHj := ^ijVfioHj, 7 := ^ pij := aij{^)+Pij. (3.7) 

/i + /io 

is defined in fj2.22|) . and {i',a-)E := {iy,aE), [z/, a-] := [z/, ai?]. 
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